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f) j' Abstract. PICsIT (Pixellated Imaging CaeSium Iodide Telescope) is the high energy detector of the IBIS telescope 

on-board the INTEGRAL satellite. PICsIT operates in the gamma-ray energy range between 175 keV and 10 
m : MeV, with a typical energy resolution of 10% at 1 MeV, and an angular resolution of 12 arcmin within a ~100 

qq ■ square degree field of view, with the possibility to locate intense point sources in the MeV region at the few arcmin 

level. PICsIT is based upon a modular array of 4096 independent CsI(Tl) pixels, ~0.70 cm 2 in cross-section and 
3 cm thick. In this work, the PICsIT on-board data handling and science operative modes are described. This 
^ ' work presents the in-flight performances in terms of background count spectra, sensitivity limit, and imaging 

capabilities. 
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The IBIS telescope (Ubertini et al. 2003) provides unprecedented imaging capabilities over three decades of energy 
Qh intervals in the gamma-ray regime, using a composite two-layer detection system made of discrete independent pixels. 
q The upper layer, ISGRI (IBIS Soft Gamma-Ray Imager), is the low energy detector (Lebrun et al. 2003). ISGRI is 
?— i ' based on CdTe pixels and covers the energy band 15 — 1000 keV. Placed 10 cm below ISGRI, PICsIT (Pixellated 
c/3 Imaging Csl Telescope), is based on CsI(Tl) crystals, and it is operative between 175 keV and 10 MeV. Both layers 
' are position sensitive and allow the reconstruction of the shadowgram cast by the tungsten coded mask placed 320 
>- cm above ISGRI. At that distance, one mask pixel, 1.12 cm wide, subtends an angle of ~12 arcmin, which defines the 
J^j ' angular resolution of the instrument, while the point source location accuracy (PSLA) for strong sources (~ 20cr) is 
. about 1 arcmin. A passive collimation system is present between ISGRI and the coded mask to shield the detectors 
03 against the cosmic diffuse background up to about 300 keV, while both layers are in anticoincidence with an active 
BGO veto system which surronds the sides and bottom of PICsIT for high-energy background rejection (Quadrini 
et al. 2003). The two-layer focal detection plane array is also coupled with the on-board IBIS calibration unit based 
on a 22 Na source with a modified veto element as tagging detector, for gain and performance monitoring (Bird et al. 
2003a). 

PICsIT (Labanti et al. 2003) has been designed to attain, with the same instrument, good sensitivity, fine imaging, 
and moderate energy resolution in the gamma-ray energy domain between 175 keV and 10 MeV. It consists of 8 totally 
independent rectangular modules arranged in a 2 x 4 array inside an aluminium mainframe. Each module contains 
512 (16 x 32) CsI(Tl) scintillation crystals, each one coupled with a silicon photodiode (PD) for fluorescence light 
collection. The total assembly is a 4096 (64 x 64) pixel matrix. Each pixel is 0.84 x 0.84 cm 2 in cross-section, and 3 cm 
thick. Every crystal is polished and wrapped with a white diffusive paper to provide stability and physical isolation 
between pixels. Each PD is interfaced to an analogue front-end electronic chain including a biasing filter, charge 
amplifier, shaping amplifier, threshold discriminator, and peak detector plus stretcher. These circuits are implemented 
in an analogue application specific integrate circuit (ASIC). Each ASIC receives the signals from 16 PDs. The key 
performance parameter of PICsIT are listed in Table 1. 

Send offprint requests to: G. Di Cocco - e-mail: dicocco@bo.iasf.cnr.it 

* Based on observations with INTEGRAL, an ESA project with instruments and science data centre funded by ESA member 
states (especially the PI countries: Denmark, France, Germany, Italy, Switzerland, Spain), Czech Republic and Poland, and 
with the participation of Russia and the USA. 



2 Di Cocco et al.: IBIS/PICsIT in-flight performances 

Table 1. IBIS/PICsIT key performance parameters. 



Parameter 


Value 




Notes 


i y v l ^ c l i w ± l ccl 
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9% at 1275 keV 
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(Revolution 38) 


650 cts/s 




multiple events 


Angular Resolution 


12' 




sampled in 10' pixels 


PSLA 


5' 






Continuum Sensitivity 


~ 2 x 10~ b photons cm -2 s~ 


1 keV" 1 


at 500 keV 


(3a, 10 5 s, AE=E/2) 


~ 1 x 10 _ photons cm - s~ 


1 keV" 1 


at 2 MeV 


Line Sensitivity 


~ 2.5 x 10~ 4 photons cm" 


1 s- 1 


at 500 keV 


(3cr, 10 5 s, AE=E/2) 


~ 2.6 x 10~ 3 photons cm" 




at 2 MeV 



The PICsIT detector has undergone a long series of functional tests and scientific calibration before and after the 
the integration in the IBIS telescope. The pixel-by-pixel channcl-to-energy relationship was calibrated on-ground in 
order to create the look-up tables (LUT) uploaded in the on-board data handling system (Labanti et al. 2003). The 
LUTs have then been re-created using the in-flight data from the IBIS calibration unit (Bird et al. 2003), and the 
first in-flight calibration results indicate a very good stability of the behaviour of PICsIT (Malaguti et al. 2003). 



2. IBIS/PICsIT modes of operation and in-flight context 

The on-board telemetry limitations do not allow a complete (position, time, energy), event-by-event, transmission 
of PICsIT data. This has implied a significant amount of on-board processing. Each PICsIT valid signal 1 is first 
corrected for gain non-uniformity by means of the on-board LUT. These corrected events are then associated with their 
respective energy deposit, incidence pixel, and time tag. This on-board event pre-processing is performed differently 
according to the number of pixels triggered, the selected PICsIT operative science mode, and depends finally on the 
status of an on-board set of parameters which can be updated and uploaded via telecommand. 

2.1. Event multiplicity and handling 

In order to minimize the impact on scientific performances caused by the strong limitation imposed by the telemetry 
budget at satellite level, PICsIT data are treated on-board in different ways depending upon their multiplicity. A 
single event is the result of a gamma-ray photon triggering only one pixel, while a multiple event triggers two (double) 
or three (triple) pixels. Events with multiplicity greater than 3 are not transmitted to ground. Singles and multiples 
are put in two independent data sets. Multiple events are then processed on-board by the data handling system (DHS) 
and associated with one single energy deposit and incidence pixel. The first is the sum of the energies deposited in 
each pixel, while the incidence pixel, to form the multiple event shadowgram to be coupled with the mask pattern for 
image deconvolution, is reconstructed on-board by the DHS, which calculates the most probable incidence position 
implementing an algorithm developed from MonteCarlo simulations and scientific calibrations (Malaguti et al. 2001). 
The on-board event handling introduces a certain degree of detector non-uniformity. This is due to the fact that 
multiple events are treated at semi-module level (16 x 16 pixels), and, since there is no electronic communication 
between semi-modules, energy deposits falling across semi-module boundaries are not recombincd into a single event. 
This gives rise to characteristic bi-convex (concave) shapes of the background for multiple (single) events (Malaguti 
et al. 2001, Bird et al. 2003). 



1 Signals are not treated as PICsIT events when they are in coincidence with ISGRI ones (Compton events), calibration unit 
tags (calibration events), or veto strobes (removed as background). 
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2.2. PICsIT operative scientific modes 

PICsIT has two fundamental operative modes: Photon-Photon Mode (PPM) and Standard Mode (SM). Because of 
the already mentioned satellite telemetry limitations, PPM can be used only in the case of limited count rate (e.g. 
selecting only one region of the operative spectral range). In PPM, the events (singles and multiples) are transmitted 
with full information: energy (1024 channels), pixel position, time tag between events. During slews, PICsIT is by 
default, switched to PPM. 

SM is the nominal operating mode for PICsIT and it is composed of two submodes: Spectral Imaging (SI) and 
Spectral Timing (ST). In SI the events (singles and multiples) are integrated on-board for a characteristic time which 
is coupled to the duration of one pointing, i.e. ~2000 s, to produce 3-dimensional histograms. SI histograms are 
optimized for spatial resolution. The raw 1024 channel spectra are rebinned to 256 channels, and one energy spectrum 
is obtained for every pixel. The final result is a 256 x 64 x 64 histogram. The spectral re-binning is performed by 
means of a non-linear binning table (contained in the DHS context) in order to optimize the energy sampling. 

On the other hand, ST is optimized for temporal resolution at the expense of the spatial information. PICsIT 
and the on-board DHS have been designed to have up to 8 programmable energy bands, and down to 0.97 ms 
time resolution. The implementation of ST in full resolution would however imply, because of telemetry limitation, 
significant loss of data for IBIS as a whole. The default (from end of June 2003) functioning foresees four energy bands 
(~ 150 200 keV, - 200 + 250 keV, - 250 350 keV, ~ 350 700 keV) with 1.9 ms time resolution. 

2.3. PICsIT and DHS in-flight context 

PICsIT in-flight data processing requires a number of parameters and algorithm to be present and implemented on - 
board. The main ones are: LUT for gain and offset correction, channel-to-energy conversion, binning tables for SM, 
incidence pixel reconstruction algorithm from multiple events, dead pixel map. These are all stored in the PICsIT and 
DHS contexts. The contexts can be downloaded as part of the scientific housekeeping data (HKD) and their contents 
can be updated and uploaded via telecommand following a "reboot" of PICsIT. 

A typical example is the handling of noisy pixels. When a pixel becomes noisy, it is killed by the on-board automatic 
noisy pixel handling system. If it remains noisy, its coordinates are put in the killed pixel map contained in the context 
and the pixel will remain switched off at the first PICsIT reboot. Only 2 pixels have become noisy during the first 7 
months of the mission, and have therefore been killed, for a total of 51 killed pixels (~1% of total detecting area). 

3. In-flight performance and sensitivity 

The level and stability of PICsIT scientific performance has been monitored during the first months of the mission 
using the data from the in-flight calibration unit. For the first 4 months of the mission, PICsIT was switched off during 
belt transits. This implied a temperature decrease of ~10°C at constant satellite attitude during belt exits. These 
temperature fluctuations did not have measurable influence on energy resolution (measured at the two 22 Na lines: 511 
and 1275 keV), which remained constant. On the other hand, the observed gain drift during temperature increase 
and stabilization at belt exits up to revolution 37 are fully explained by temperature fluctuations (Malaguti et al. 
2003). Therefore, after having verified the absence of significant effect on instrument behaviour, leaving few modules 
on during belt transits for several weeks, it was decided, starting from revolution 38, to keep PICsIT always switched 
on to avoid unnecessary temperature fluctuations and increase the instrument stability. The observed residual gain 
fluctuations are to be ascribed to the changing aspect of the satellite with respect to the sun direction, for the various 
pointings. 

3.1. PICsIT background count spectrum 

Prior to IBIS/PICsIT, no other pixellated Csl detector of this size was flown on a satellite. For this reason, the in-orbit 
background level and energy spectrum predictions, based on MonteCarlo simulations, or on semi-empirical scaling 
from other high energy space telescopes (e.g. OSSE, SIGMA), were accurate only within a factor of ~2. 

A series of photon-particle transport simulation programs have been developed to evaluate the in-orbit background 
rate and spectral shape. These resulted in a predicted total (all events) count rate between 1500 and 2000 counts/s 
(Armstrong et al. 1999, Lei et al. 1997). However, PICsIT design assumed a conservative estimate of 7500 counts/s at 
solar maximum level in the energy range 140 keV - 10 MeV. The most important contigency included in this evaluation 
is the uncertainty in the incident cosmic ray proton flux (Ferguson et al. 2003). Figure [T] shows the PICsIT in-orbit 
background spectrum for single and multiple events compared with the spectral shape obtained from Monte-Carlo 
simulations and scaled up to 7500 counts/s. The measurements are from the empty field acquisitions performed during 
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Fig. 1. PICsIT in-orbit background spectrum for single and multiple events compared with the expected expected 
spectrum obtained from MonteCarlo simulations normalized to 7500 counts/s. 



Table 2. IBIS/PICsIT parameters for the calculation of the statistical continuum flux sensitivity limit. 



Parameter 


Description 


Units 


Notes 


pUont 
-^min 


minimum detectable continuum flux 


photons cm -2 s _1 keV 


-l 


mm 


minimum detectable line flux 


photons cm~ z s _1 






statistical significance 




3 


£ 


total detection efficiency 




MonteCarlo and ground calibrations 


EP 


photopeak detection efficiency 




MonteCarlo and ground calibrations 


T 


observing time 


s 


10 b s 


AE 


Energy bin 


keV 


see Figure |21 


SE 


Energy resolution 


keV 


19% at 511 keV, 10% at 1 MeV 


e c 


Closed mask element transparency 




~6% at 400 keV, ~31% at 3.5 MeV 


e D 


Open mask element transparency 




1 


i 


Coding power 




0.73 (a) 


eipr 


Incident pixel reconstruction efficiency 






B 


Background count rate 


counts cm -2 s _1 keV" 


1 in-flight data 



(a, In the case of position sensitive detectors coupled with coded masks, mask and detector have in general different spatial 
sampling, and the coding power I is defined by: I — 1 — (Ld/3-Lm) where Ld and Lm are the detector and mask pixel 
dimensions respectively (Skinner 1995). For PICsIT we have Lu/Lm = 11.2/9.2, and therefore / = 0.73. 



revolution 38 (just before the first Crab observation), and give a total counting rate of ^2900 and ~650 counts/s for 
single and multiple events, respectively. 

It is worth noting that, while the total count rate is a factor of ~2 less than the expected conservative value, the 
observed background flux up to ~250-300 keV is higher than expected. This is partly due to the events generated by 
the high-energy tracks to be ascribed to the photon-particle shower generated by high-energy cosmic rays interacting 
in the satellite-telescope structure or within PICsIT itself (Segreto et al. 2003). 
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Fig. 2. PICsIT minimum detectable flux for continuum (left panel) and narrow line (right panel) emission as a function 
of energy assuming only statistical effects, and the in-flight measured background. The dashed lines of the first two 
bins in the left panel indicate the worsening due to the effect of cosmic-ray induced tracks. 



3.2. Continuum and line sensitivity 

The statistical limits, i.e. the minimum detectable fluxes in case of negligible systematic effects, of PICsIT sensitivity 
to continuum (F^™* in photons cm~ 2 s _1 keV -1 ), and narrow 2 line (F^!™ in photons cm -2 s _1 ) emission, are given 
by: 

= £ . A . T . A /' ( e1-e c )2. gipr ( nAe ° + 9c) + ■ (e ° + 0c)2 + 4 ■ (e ° + 9c)2 ■ A ■ T ■ AE ' B ) (1) 
= j-= j-^ik — ( n ° (©o + e c ) + Vnl ■ (e + e c ) 2 + 4 • (e + e c ) 2 -a-t-se-b) (2) 

0.75 • ep ■ A ■ T ■ (Bo - B c ) 2 • £i pr V J 

Since, both in PPM and SM, PICsIT single and multiple events arc contained in two different and independent 
data sets, the total statistical limit for the minimum detectable flux is given by: 

-^min = I 1 l (3) 

W ^singles -j 2 ^multiplcs^2 

The input parameters for equations ^ are described in Table 13.11 The background spectrum is shown in figure ^ 
and discussed in section ETT1 the transparency of the closed elements of the coded mask has been calculated from the 
cross-section of 16 mm of tungsten, while the photopeak and total detection efficiency for single and double events 
has been obtained by means of a MonteCarlo photon-particle transport code cross-checked with the IBIS on-ground 
calibration results (Bird et al. 2003b). 

Figure |21 shows the continuum and line PICsIT sensitivity curves, based on statistical principles only and calculated 
using the in-flight measured background spectrum. An important caveat is that the actual sensitivity limit depends 
also on systematic effects, mainly due to non-uniformity of the spatial response for both background and source events 
(Natalucci et al. 2003), to cosmic-ray induced events (Segreto et al. 2003), and to background temperature variations 

2 In this context an emission line is considered narrow, when its width is small compared with the energy resolution of the 
instrument at that energy. 
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Table 3. Actual source position accuracy. Columns: (1) Source name; (2) configuration during the observation and 
offset with respect to the centre of the field of view; (3) Reference position (RA, Dec, J2000) from catalog or circulars; 
(4) PICsIT position (RA, Dec, J2000); (5) Offset [arcmin]. The position for PICsIT are taken in the 252 - 329 keV 
energy band, except for GRB021125 which has been obtained for 203 — 252 keV, where the significance of the detection 
was the highest. 



Source 


Configuration Reference Position 


PICsIT Position 




Offset 


(1) 


(2) (3) 


(4) 






(5) 




Plain deconvolution (brightest pixel as 


reference) 








Crab 


staring, on-axis 05 : 34 : 32; +22 : 00 : 52 


05 : 34 : 13; +21 


55 


48 


6.8' 


GRB 021125 


dithering, 7.3° off-axis 19 : 47 : 56, +28 : 23 : 28 


19 : 47 : 34, +28 


15 


13 


9.5' 


Calculation of the centroid 


Cygnus X-l 


dithering, on-axis 19 : 58 : 22; +35 : 12 : 06 


19 : 58 : 08; +35 


08 


40 


4.4' 


Crab 


dithering, on-axis 05 : 34 : 32; +22 : 00 : 52 


05 : 34 : 24; +22 


01 


20 


1.9' 


Crab 


dithering, 30' off-axis 


05 : 34 : 17; +21 


59 


11 


3.8' 


GRB 021125 


dithering, 7.3° off-axis 19 : 47 : 56, +28 : 23 : 28 


19 : 47 : 51, +28 


19 


16 


4.3' 



Table 4. PICsIT observations of the Crab during Rev. 39 (staring, on axis, 77 ks), Rev. 43 (dithering, on axis, 132 ks), 
Rev. 44 (dithering, off axis 0.5°, 104 ks). The columns indicate the energy band [keV], the count rate in the peak pixel 
[counts/s] and the significance of the detection for each revolution [a]. 







Rev 39 




Rev 43 




Rev 44 


Energy Band 


Rate 


Significance 


Rate 


Significance 


Rate 


Significance 


200 - 252 


2.5 


7.4 


1.5 


5.9 


3.1 


5.7 


252 - 329 


2.2 


11.0 


2.1 


10.0 


2.6 


8.0 


329 - 455 


1.3 


7.1 


1.2 


6.5 


1.4 


5.0 


455 - 655 


0.7 


4.4 


0.6 


3.3 


0.6 


3.0 



(Stephen et al. 2003). The actual version of the INTEGRAL Data Analysis System (IDAS 3 ) performs the background 
correction on the basis of an empty field observation (Goldwurm et al. 2003). A refined correction of these effects will 
be possible with a background model obtained after an accurate study of the in-flight data. 

3.3. Imaging capabilities 

During the Performance Verification and Calibration Phases, IBIS has observed two point sources: Cyg X-l and the 
Crab nebula. The first results shown here have been obtained with the Offline Scientific Analysis (OSA) of the IDAS, 
and show the imaging capabilities of PICsIT. Figure|3and0]show the significance map and profile in four energy bands 
(only one for Cyg X-l), obtained decoding the shadowgram casted by the coded mask (see Goldwurm et al. 2003 for 
details regarding image deconvolution in IBIS and the definition of significance). The significance of a detection in 
a coded mask telescope depends on a number of factors like, among others, the ratio between detector plane spatial 
resolution and mask element size, background non-uniformities, and off-axis position of the source (Skinner 2002). 
In the current version of IDAS, the calculation of the significance takes into account statistical effects plus a first 
evaluation of background and detector non-uniformity (using the empty field observation of revolution 38), while an 
interpolation is implemented to account for the uneven sampling of mask elements by the detector pixels. 

Table 3 shows the reconstructed sky coordinates of the three point sources observed with PICsIT, compared with 
the reference position, which has been obtained from catalogues (Crab and Cyg X-l) and from ISGRI (GRB021125). 
Since PICsIT operates in a energy band where the signal-to-noise ratio (SNR) is intrinsically low, very long exposures 
are required. Therefore, it is not possible to make precise astrometry on a statistical basis. The residual offset can give 
an estimate of PICsIT point source location accuracy (PSLA), being the nominal value of the error radius of 5'. A 
more detailed evaluation of PSLA (see Gros et al 2003 for more details) will be performed with the availability of a 
refined point spread function (PSF). 

Table 4 shows the reconstructed count rates and associated significance of the Crab observations in three different 
telescope configurations. In revolution 39, IBIS observed the Crab on-axis, staring on the source for 77 ks; in revolutions 
43 and 44 the Crab was pointed with a 5 x 5 rectangular dither pattern (Jensen et al. 2003), on-axis for 132 ks and 
30 arcmin off-axis for 104 ks, respectively. 



3 IDAS software is available from http://isdc.umge. ch/index.cgi?Soft+dowload . See also Goldwurm et al. (2003) for details. 
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Fig. 4. Cyg X-l significance image and profile in the 252 — 329 keV band obtained from a 80 ks observation in 
revolution 11. 



The high fluctuation of the detected Crab rate in the lowest energy band is due to the presence of cosmic ray 
induced events. In fact, in SM the on-board accumulation of histograms allows only a partial a posteriori rejection of 
these events (see Goldwurm et al. 2003 for details on this type of correction). 

Table 4 also shows that, for the 252-329 keV band, an increase of exposure does not result in the expected increase 
of detection significance (11a with 73 ks staring, against IOct with 132 ks dithering). This is explained by the fact that 
during revolution 43, the Crab was observed in dithering mode with the rectangular 5x5 pattern (Jensen et al. 2003). 
This pattern samples regions close to the edge of the fully coded field of view, thus implying a sensitivity degradation 
of -20% (Foschini 2003). 

These results are in good agreement with simulation results obtained by means of a dedicated MonteCarlo code 
in the energy range 252-655 keV, while are a factor of two lower than expected in the first energy channel, where the 
systematic effects caused by the cosmic-rays induced events affects the data analysis. 



3.4. Spectral & Timing capabilities 

The spectral and timing characteristics of PICsIT were analysed by using other sources. A solar flare, occurred on 
November 10 th , 2002, at about 10 : 00 UTC, allowing the detection of several emission lines: the e+e~ annihilation 
line at 511 keV with the tail at lower energies made by the fusion lines, the neutron capture at 2.2 McV, and the 
nuclear excitation lines of 12 C at 4.4 McV. These lines have confirmed the linearity of the detector over the entire 
energy range, while a more detailed analysis will be performed when the refined redistribution matrix (RMF) will be 
available. 

The GRB of November 25 th , 2002, allowed a full test of the imaging, spectral, and timing capabilities in photon- 
by-photon mode and it is described in Malaguti et al. (2003). 
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In addition, the absolute timing performances both in standard and photon-by-photon mode were analysed by 
using the Crab pulsar observations (see Kuiper et al. 2003). 

4. Conclusions 

The in-flight functional and scientific performances of the IBIS/PICsIT detector have been described and discussed, 
together with the first observations performed during the PV and calibration phases. The instrument has demonstrated 
to work nominally, from an engineering point of view with respect to the design requirements. Activities are ongoing 
to improve the current scientific performances by minimizing the systematic effects due to the background evaluation, 
cosmic-rays induced events, and detector non-uniformities. 
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